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1. Introduction 

The soft part of the multi-particle production in pp collisions at high energy, which can not be 
directly described in terms of the perturbative QCD, is usually described in the framework of string 
models [ffl]-[H], which originate from the Gribov-Regge approach. In these models at first stage the 
colour strings are stretched between the projectile and target partons. The formation of the pair of 
strings corresponds to the cut pomeron in the Gribov-Regge approach. The hadronization of these 
strings produces the observed hadrons. 

In the case of nuclear collisions, the number of strings grows with the growing energy and the 
number of nucleons of colliding nuclei, and one has to take into account the interaction between 
strings in the form of their fusion and/or percolation [§]-[0]. The fusion process results in the 
reduction of total multiplicity of charged particles and growth of transverse momentum, that was 
confirmed later [|| ^] in comparison with RHIC data. 

The possible experimental observation of the string fusion or percolation is extremely inter- 
esting. Therefore, the investigation of long-range correlations (LRC) between multiplicities and 
transverse momenta in two separated rapidity intervals were proposed as the main tool to study this 
phenomenon []l0|]-[|l7|]. 

Up to the present the string fusion model was usually used for describing the multiplicity and 
mean transverse momentum of charged particles and their correlation only in the case of nucleus- 
nucleus collisions. However, the experimental data indicates an increase with energy of the mean 



transverse momentum and its correlation with multiplicity also in pp collisions [|18|]-[]21[]. In 
present work we formulate the simple model which enables to take into account the effect of colour 
string fusion on the LRC between multiplicities and transverse momenta in pp interactions. 

2. Formulation of the model 

To take into account the effect of string fusion in pp collisions one needs to know the distri- 
bution of strings in the transverse plane at given value of the impact parameter b. We'll do this in 
an analogy with the case of nucleus-nucleus collisions. 

2.1 Distribution of strings in the transverse plane: AA interactions 

In the case of AA interactions one usually assumes that at high energy a number of primaiy 
formed quark-gluon strings is proportional to the number of binary inelastic nucleon-nucleon in- 
teractions in a given event of nucleus-nucleus scattering |2|]. In frame of the classical Glauber 
model mean number of inelastic AW interactions in AB scattering at a fixed impact parameter b is 



given by the expression (e.g. see [ |22[ [23[]): 

(N, r (b)) ~ (NZ n (b)) = AB—^j— / T A {s- b/2)T B (s + b/2)d 2 s. (2.1) 

Here c^y - the cross-section of inelastic nucleon-nucleon interaction, o AB {b) - the probability of 
interaction of two nuclei at a given impact parameter b with at least one inelastic AW interaction. 
An integral of o AB {b): 

Oab = I o A B(b) d 2 b (2.2) 
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gives a so-called «production cross section» of AB interaction. The (...) in equation (2.1 ) means an 



averaging over all events with a given b. T A and 7g - the profile functions of the colliding nuclei: 



+00 



T A (s) = J p A (s,z)dz, (2.3) 

— 00 

where p A (j) - the nucleon density of nucleus A, normalized to unity, r = (s,z). The ?is a two- 
dimensional vector in the impact parameter plane. 



By (2.1) in the nucleus-nucleus collision at the impact parameter b for the string density in 



transverse plane at a point ?we have: 

o' n 

d{N str {b))/d 2 s = w str (s,b) ~ w'? oll (s,b) =AB^-T A (s-b/2)T B (s + b/2) (2.4) 

Note that in the transverse plane we place the origin in the middle between the centers of colliding 
nuclei. 

2.2 Distribution of strings in the transverse plane: pp interactions 



By (2.4) it is natural to suppose that in the case of proton-proton collision at the impact pa- 
rameter b the string density in transverse plane at a point s is proportional to 

Wst r (s,b) ~ T(s-b/2)T(s + b/2)/o pp (b) , (2.5) 

where now the T(s) is the partonic profile function of nucleon. The o pp {b) is the probability of 
non-diffractive pp interaction (with at least one cut pomeron) at a given impact parameter b. An 
integral of o pp (b): 

°pp = J Vpp^)^ (2-6) 

gives a non-diffractive cross-section of pp interaction. 

By analogy with light nuclei we will use for the partonic profile function of nucleon the sim- 
plest gaussian distribution: 

Substituting p?7| ) in ( |2~3| ) one gets 

w str (s,b) ~ e-^W/Je-P-iW/*/^,,) = e- 2 ° 2 / a2 e- h2 / 2a2 /o pp (b) . (2.8) 

We see that in this approximation the dependencies on b and s are factorized and after integration 
on ?one gets for the mean number of strings in the pp collision at the impact parameter b: 

{N str {b))~e- b2 ' 2a2 /o pp (b). (2.9) 

Since in this approach the formation of each pair of strings corresponds to one cut pomeron, 



Nstr = 2N, where ,/V is the number of cut pomerons in a given event, hence (2.9) leads to 



(N(b))~e- b2 / 2a2 /G pp (b), (2.10) 

which gives the dependence of the average number of cut pomerons on the impact parameter in 
non-diffractive pp collisions. 
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2.3 Event-by-event fluctuations of the number of cut pomerons at given impact parameter 

For the calculation of LRC one should know not only the mean number of pomerons in pp 
collisions at a given impact parameter b, but also the event by event distribution of the number of 
pomerons around this mean value. We assume that this distribution P(N,b) at a given value of the 
impact parameter b atN > 1 : 



P(N,b)=P(N,b)/[l-P(0,b)] 
is the simple modification of the poissonian distribution: 



P(N,b) 



-N(b) 



N(bf/N\ 



(2.11) 



(2.12) 



with some parameter N(b). The difference of our distribution P(N,b) ( |2.11| ) from the poissonian 
one ( 2.12 ) is only in excluding of the point N = from it: P(0,b) = 0, which corresponds to the 
absence of the non-diffractive scattering at N = 0. Clear that at N > I this reduces only to the 
introduction of the additional common normalization factor in ( [2.1 1[ ), which enables the proper 
normalization: Ljv=i = 1- 

The calculation of the mean number of pomerons at a given b with the distribution ( 2.11 ) gives: 

(N(b))=N(b)/[l-P(0,b)]- (2-13) 

Since the probability o pp (b) of the non-diffractive pp interaction at the given fixed impact 
parameter b is equal to the probability to have a nonzero number of cut pomerons, then 



o p Jb) = 1-P(0,b) = 1 -exp(-A^)) 



(2.14) 



Comparing now the formulae ( [2. 10 ) and ( |2.13 ) with taking into account ( 2.14 ) we see that in our 
model N(b) introducing a parameter Nq: 



(2.15) 



(2.16) 



N(b)=N e- h2 / 2a2 

Then the mean number of cut pomerons at an impact parameter b is given by 

(N(b))=N(b)/[\-exp(-N(b))} , 
where the N(b) is given by ( |2.15| ). 

2.4 Integration over impact parameter - min.bias pp collisions 

It is convenient to introduce in the impact parameter plane a density of the probability of 
non-diffractive pp interaction (see {2J))) normalized to unity: 



f(b) = o pp (b)/cj t 



pp 



f(b)d 2 b=l. 



(2.17) 



Using the formula (2.17) one can find a mean number of pomerons in non-diffractive pp interaction 
averaged over the impact parameter: 



(N) = / (N(b))f(b)d 2 b = / N(b)d 2 b/a pp = 2Ka 2 No/a, 



pp 



(2.18) 
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and the corresponding variance D N = {N 2 ) - (AO 2 , where 



{N ) = j (N 2 (b))f(b)d 2 b = xa 2 N (N + 2)/o pp = (N) (N + 2)/2 (2.19) 

Here we have used ( 2.15| ). In our model using ( [2.1 5[ ) we can calculate also the cross-section of 
non-diffractive pp interaction 



pp 



o p Jb)d 2 b = / [l-ex V {-N{b))]d 2 b = 2Ka 2 <S> { {N ) , 



where 



<Mz)= / 



dt 



(_1)*+V 



k=\ 



klk" 



(2.20) 
(2.21) 



2.5 Probability of A^ cut pomerons in non-diffractive pp collision 



In the framework of our simple model we can find the probability w N to have N cut pomerons 
in a non-diffractive pp collision by averaging the P(N,b) ( 2.11| ) over b at fixed N: 

P(N,b)f(b)d 2 b = — 



PP 



P(N,b)d 2 b , 



where we have taken into account ( |2.14] ). Using now ( 2.12 ), we have 



w N 



a pp N\ 



- N ^(N(b)fd 2 b 



2n 



- N{b \N{b)f bdb 



One can introduce in ( 2.23[ ) the new integration variable N in accordance with (2.15): 

N = N{b) = N e- h2/2al , dN = -(N/a 2 ) bdb . 



Then (2.23) takes the following form 

w N = 



o pp N\ 



Such integral is a difference of the gamma and incomplete gamma functions: 



A'o 



e-^~ l dz= e-^dz- e- z i*- l dz = T{N)-T{N,N Q ) , 



At integer N 



N-l 



F(N) = (N-l)\, r(N,N ) = (N— l)le~ N ° £ N l Jl\ 

1=0 



Gathering we find 



w N 



2%a 2 



1=0 



PP 



where we have introduced the o N by 



2na z 



N 



l-e^N'/l 

1=0 



The direct summing gives 



£ a N = 2xa 2 ® x (No) = o 1 



N=l 



PP 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



where we have used ( 2.20 ) and ( |2.21[ ). Recall that o pp is the non-diffractive pp cross section (2.6). 
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2.6 Comparison with quasi-eikonal and Regge approaches 



Now we see that our formula for the o N (2.29) coincides with the well known result for the 
cross-section o N of N cut-pomeron exchange, obtained in the quasi-eikonal and Regge approaches 

4nX 



[24,25, 



CN 



N-l 



1 



I z k /kl 



k=0 



where 



^exp(A^), X=R 2 + a'$, $ =\n{s/\GeV 2 ) . 

A 



(2.31) 



(2.32) 



Here A and a' are the residue and the slope of the pomeron trajectory. The parameters y and R 
characterize the coupling of the pomeron trajectory with the initial hadrons. The quasi-eikonal 
parameter C is related to the small-mass diffraction dissociation of incoming hadrons. 

This enables to connect the parameters No and a of our model, which describe the dependence 
of the mean number of pomerons on the impact parameter b (see formulae (2.15) and (2.16)) with 
the parameters of the pomeron trajectory and its couplings to hadrons. Comparing ( 2.29 ) and ( 2.3 1 ) 
we have 



M> = z = ^exp(A£), « = y^, X=R 2 + a'^ 



(2.33) 



In our calculations the numerical values of the parameters for the case of pp collisions were taken 
from the paper [ 27 ] : 

A = 0.139, a' = 0.21 GeV~ 2 , 7=1.77 GeV~ 2 , R 2 = 3.18 GeV~ 2 , C = 1.5, (2.34) 



which gives, for example, for the parameters A^o and a entering in our formula ( |2. 15[ ): 

a = 0.51fm, Af = 3.38 at^ = 60GeV, 
a = 0.60fm, No = 9.02 at sfs = lOOOGeV . 



(2.35) 



3. String fusion effects 



As follows from the consideration in the previous section, in pp interactions the number of 



contributing cut pomerons (see formulae ( J2. 15D and ( |2.16[ )) increases with the collision centrality 



and energy ( |2.35| ). In our approach the number of strings is twice the number of cut pomerons, so it 
also will increases with the collision centrality and energy. Since the strings have a certain limited 
size in the transverse plane (a plane of the impact parameter) an overlap of the strings will start 
with increase of their density. As a result the color fields of different strings will interact what will 
influence on their fragmentation process. For taking into account the effects from the interaction 
of strings at high density a string fusion model (SFM) was proposed 



3.1 Different versions of string fusion model 



There are two versions of the SFM a model with local fusion (overlaps model) [28] and a 
model with formation of global color clusters (clusters model) [E9h. In both versions the fusion 
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SFM 



cellular 
analog 
of SFM 



"overlaps" 
(local fusion) 



o 



QQ 



c={s,,s,,...} 

Sh - area 
covered k-times 



00 




□ 



c={n*;} 











N 

















I _ \jsrr _ "occupation" 
// ~~ ij numbers 



"clusters" 
(global fusion) 




S:=3a ; N;=5; k =5/3 
S=2a o ; n;=4; k^=2 



Figure 1: The different versions of the SFM. The first column corresponds to the version with a local 
string fusion (overlaps model). The second column corresponds to the formation of fused string clusters 
(clusters model). The first row corresponds to the original versions of the SFM [ |28[ p9| | and the second row 
corresponds to their cellular analogs [|o] |T[ |5] Q| (see text for details). 



processes result in the reduction of total multiplicity of charged particles and growth of transverse 
momentum. 

In the first valiant according to [^8|] we assume that the mean multiplicity per unit of rapidity 
and the mean transverse momentum of charged particles emitted from the region, where k strings 
are overlapping, are described by the following expressions: 



(n) k = HoVkS k /(jQ 



Po 



2 Vk, 



1,2,3,.. 



(3.1) 



nrj tr is the 



Here St is the transverse area of the region, where k color strings are overlapped. Ob 
transverse area of a string. }1q and p are the mean multiplicity per unit of rapidity and the mean 
transverse momentum of charged particles, produced from a decay of one string. 



In the second variant one assumes that the fused strings form a cluster. According to [29] 
in this case we suppose that the mean multiplicity per unit of rapidity and the mean transverse 
momentum of charged particles emitted by a cluster with a transverse area S c i, formed by k strings, 
can be found as follows 



(n) c i = lM)V k ciSci/o , 



(P?)d 



2 



koo/S c i 



(3.2) 



Note that in two limiting cases (an absence of overlaps of strings and a total overlapping of string 
areas) both variants lead to the same rules. 

Later the simplified discrete analogs of both mentioned versions of SFM based on the imple- 



mentation of the lattice in the transverse plane [30, 31, 32] were proposed. It was demonstrated 
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[ 33 , [34]] that the MC algorithms of the calculations of LRC coefficients based on the cellular ver- 
sions of SFM work much faster and give practically the same results as the ones based on the 
original versions of SFM. These cellular versions of SFM enable also to calculate analytically the 



mean values of the observables and LRC coefficients in some limiting cases [ pl[ , |34[ ], what was 
used to control the reliability of the created MC codes. 

In Fig.[j] we illustrate visually the prescriptions of all these versions of the SFM. The first 
column corresponds to the version with a local string fusion (overlaps model). The second column 
corresponds to the formation of fused string clusters (clusters model). The first row corresponds to 
the original versions of the SFM and the second row corresponds to their cellular analogs. 

3.2 Monte-Carlo algorithm for the calculation of LRC functions 

The calculation of LRC functions (regressions) is based on the formula 

Ec>v(C) (n B ) c P c {n F ) 



L c w(C) P c (n F ) 



(3.3) 



obtained in [ |13[ ] for the correlation between multiplicities n F and n B in separated rapidity windows 
and on the similar formula 

, v _ L C w{C) (p tB )c p c{n F ) n „ 

\PtB 7n„ ~ ^ 77a d 7 \ ^ ' 

Lc^(C) P c {n F ) 

used in [ ^Ij |32| , 33] for the calculations of the LRC function between the multiplicity n F in the 
forward rapidity window and the corresponding mean transverse momentum p tB of n B charged 
particles in the backward window. 

The calculations of the sums on string configurations C in numerator and denominator of the 



formulae ( p.3| ) and ( p.4| ) were performed by MC simulations of the configurations C with proper 
weights w(C): 

J>(C)... = (3.5) 

C Uslm sim 

With this purpose at the generation of the string configurations C the results of the Section 2 were 
used. 

For the generated string configuration C the mean values of the multiplicity {n B ) c and trans- 



verse momentum (p tB ) c in the backward rapidity window, entering the formulae (3.3) and (3.4), 



were calculated using the prescriptions (3.1) and (3.2) for the cases of local string fusion and cluster 



formation, both in frameworks of the original version of SFM and its cellular analog. 

We also have supposed that the multiplicity distribution of charged particles produced from 
the decay of any string is poissonian, what leads to the poissonian distribution for the probability 
of production of n F charged particles Pc(n F ) from the given string configuration C: 

((n F ) c T p 



Pc{n F ) = P(„ F ) c (n F ) = e 



-<"f)c 



(3.6) 



where one can calculate (n F ) c similar to (n B ) c using the prescriptions (3.1) and (3.2) 



Similarly to the correlation functions (3.3) and (3.4) one can also calculate (by MC simulations 
of the string configurations C) the overall mean values (averaged over all events) in the forward 



(%) = L^(C)(?v)c , (Pif) = J^w{C)(PtF)c 
c c 



(3.7) 
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Figure 2: The long-range correlation (LRC) function in pp collisions between multiplicities n F and n B in 
separated rapidity windows calculated by (3.3) at energies 63 and 900 GeV for the forward and backward 
rapidity windows of the same width Ay F = Ay B = dy = 0.8, what corresponds to the mean multiplicities 
(n F ) = (n B ) = (n)=1.5 and 3.0 per window at these energies. Points O and • - the results of calculations in 
the framework of the original versions of the string fusion model (SFM) with the local fusion (overlaps) or 
the formation of string clusters (clusters) correspondingly (see first row in Fig.^J) p8[ [l9| |. Points □ and ■ 
- the results of calculations in the framework of the lattice analog of SFM with the local string fusion or the 
cluster formation correspondingly (see second row in Fig.|J) [ |30| , pT[ f3~3[ [34]]. 



and (by the same way) in the backward windows. 

Note that the resulting distribution on the number of charged particles produced in the forward 
or backward window 

P(n F )=£w(C)P c (n F ) (3.8) 



will be non-poissonian despite the poissonian form (3.6) of the distribution Pc{n F ) and almost 



poissonian (see (2.1 1) and (2.12)) fluctuations in the number of cut pomerons at fixed value of the 
impact parameter b in our model. The reason is the non-poissonian fluctuations in the number of 
cut pomerons originating from the event-by-event fluctuation of the impact parameter b (see also 



4. Results 

In Figs. |2| and || we present as an example the results of our calculations of the LRC functions 
in pp interactions using the MC algorithm based on the model described above. 

In Fig|2] the results of calculations of LRC between multiplicities n F and n B in separated ra- 



pidity windows on the base of the formula (3.3) at the initial energies 63 and 900 GeV are pre- 
sented. In both cases the rapidity width of the forward and backward windows is taken the same 
Ay F = Ay B = dy = 0.8, what corresponds to the mean multiplicities (n F ) = (n B ) = (n)=1.5 and 3.0 
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Figure 3: The long-range correlation (LRC) function in pp collisions between the multiplicity n F in the 
forward rapidity window and the corresponding mean transverse momentum p tB of n B charged particles in 



the backward window calculated on the base of the formula (3.4) at the energy 900 GeV for the forward 
and backward rapidity windows of the same width Ay F = Ay B = dy = 1 .6, what corresponds to the mean 
multiplicity (n F ) = (n B ) = («)=6.0 in the window. Notations of the points are the same as in Fig.[| Points 
x - the experimental data of the ALICE collaboration [|7]| on the correlation between the multiplicity n and 
the mean transverse momentum p t of charged particles in the same pseudorapidity interval rj € (—0.8,0.8), 
Arj = 1.6 obtained in pp collisions at 900 GeV. 



in such window at these energies. We see in Fig.^]the considerable increase of the strength of the 
LRC between multiplicities in separated rapidity windows. 

In Fig.|3] the results on the LRC between the multiplicity n F in the forward rapidity window and 
the corresponding mean transverse momentum p tB of n B charged particles in the backward window 



on the base of the formula (3.4) at the initial energy 900 GeV are presented. The rapidity width 
of the forward and backward windows is taken the same Ay F = Ay B = Ay = 1.6, what corresponds 
to the mean multiplicities (n F ) = (n B ) = (n)=6.0 in such window at this energy. In Fig.|| we also 
present the experimental data of the ALICE collaboration [|7|] on the correlation between the mul- 
tiplicity n and the mean transverse momentum p t of charged particles in the same pseudorapidity 
window Tj G (—0.8,0.8), At] = 1.6 obtained in pp collisions at the energy 900 GeV. The difference 
between calculated LRC function and the ALICE experimental data at large multiplicities can be 
explained by the contributions of additional short-range mechanisms in the case of the correlation 
between the multiplicity n and the mean transverse momentum p t in the same window, which do 
not contribute in the case of LRC. 

We see also in Figs. ^| and |3] that the results of calculations of the LRC strength in the frame- 
works of all four versions of SFM (local fusion or cluster formation and their lattice analogs, see 
Figs.[j]) turn out to be very close to each other. On the one hand, because of the principal techni- 
cal differences in the realization of these SFM versions, it indicates in favour of the reliability of 
the obtained results. On the other hand, it obviously means that one can not distinguish between 
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versions of the SFM mechanisms on the base of this comparison of the obtained results with the 
experimental data and the calculations of other string fusion effects on the physical observables are 
needed. 



5. Conclusions 

The simple model which enables to take into account the effect of colour string fusion in pp 
interactions is suggested. At the formulation of the model we assume that the dependence of the 
average number of cut pomerons on the impact parameter in a non-diffractive (ND) pp collision 



is gaussian (2.15) with the additional condition that we have at least one cut pomeron (2.16). We 
assume also that the event-by-event distribution of the number of cut pomerons around this average 



value at fixed impact parameter b is poissonian (2.12) with the same condition (2.11 ). 

It is shown that these two simple assumptions after integration over impact parameter lead 
to the well known formula for the cross-section On with N cut-pomeron exchange in a ND pp 



collision ( |2.29| ), ( |2.31| ), which was obtained in the quasi-eikonal and Regge approaches [ ]24| , \2.5[ [2q ]. 
This have enabled us to connect the parameters of our model with the parameters of the pomeron 



trajectory and its couplings to hadrons ( |2.33| ). 

The effects of the string fusion [|5| |6|] on the multiparticle production were taken into account 
in the same way as it was done in the case of AA collisions [32, ^3Q . At that the different version of 
the string fusion mechanism (local fusion or cluster formation and their lattice analogs, see Fig.[l]) 
were considered. 

On the base of the model the Monte-Carlo algorithm was developed and the long-range cor- 
relation functions between multiplicities and between the average transverse momentum and the 
multiplicity in pp collisions at different energies were calculated. It was found that the results of 
calculations of the long-range correlation (LRC) strength in the frameworks of all four versions of 



the string fusion model (SFM) turn out to be very close to each other (see Figs. 1331 and |3.4| ). 
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